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Discretized kinetic theory on scale-free networks 
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The network of interpersonal connections is one of the possible heterogeneous fac¬ 
tors which affect the income distribution emerging from micro-to-macro economic 
models. In this paper we equip our model discussed in ELS with a network struc¬ 
ture. The model is based on a system of n differential equations of the kinetic 
discretized-Boltzmann kind. The network structure is incorporated in a probabilis¬ 
tic way, through the introduction of a link density P(a ) and of correlation coef¬ 
ficients P(f3\a), which give the conditioned probability that an individual with a 
links is connected to one with (3 links. We study the properties of the equations 
and give analytical results concerning the existence, normalization and positivity of 
the solutions. For a fixed network with P(a) = c/a q , we investigate numerically 
the dependence of the detailed and marginal equilibrium distributions on the initial 
conditions and on the exponent q. Our results are compatible with those obtained 
from the Bouchaud-Mezard model and from agent-based simulations, and provide 
additional information about the dependence of the individual income on the level 
of connectivity. 


I. INTRODUCTION 


The network of interpersonal connections is an important factor in defining economic 
interactions of individuals in a society. We may suppose that well-connected individuals 
have more frequent and numerous interactions, with a positive effect on their income. Con¬ 
versely, one might argue that rich individuals can afford more connections than poor ones, 
because maintaining a link requires some expenses and investment for communication and 
transportation. We may therefore wonder if there is a relation between the link density 
distribution P(a) and the income distribution, in the sense that, for instance, to a larger 
number a of links corresponds a larger income. The assortative or disassortative nature of 
the network may also have an effect on income distribution. Common sense suggests that 
in a developed society the correlation is more of a disassortative type, namely such that 
individuals with many links are mainly connected with individuals with few links. This 
could be for instance the case of the economic relationship between the owner of a firm and 
her employees, or between the owner of a store and his clients. 

Most existing studies of economic interactions on a network do not concern interactions 
between individuals, but between financial institutions, companies or countries [iij]. When 
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networks of individuals have been considered, they have generally been treated in the context 
of interacting-agents models. In these models simulation algorithms pick individuals at 
random from an ensemble and let them interact according to certain rules. It is relatively 
straightforward to improve the algorithms by inserting into the list of agents some pre-dehned 


connections which form a certain network 


4-6]. It is even possible to let the network evolve 


dynamically, with links changing in dependence on the wealth of the nodes 


7 


111 - 


Another approach is based on multi-agent linear stochastic equations [12| [14|. These have 
the form Wi(t) = rji{t)wi(t ) + Jij w i(t), where -uy(f) is the wealth of the 

i-th agent, r]i(t) is a stochastic noise and J an interaction matrix. In this case, the adjacency 
matrix of a network can be included in J. The linear nature of the equations allows, even 
for a large number of agents, an effective analytical or numerical treatment. The drawback 
of the linear equations, however, is that they cannot completely describe the interaction 
process, since the variation of the wealth of the i-th agent depends on the wealth of the 
others, but not on his own wealth (apart from the stochastic term). 

In a discretized kinetic theory, the state variable of a system (energy, income, ...) is 
discretized into a finite number of levels or classes. The Boltzmann equation, which is based 
on the concepts of transition probabilities and detailed balance and has in general the form 
of an integro-differential pde, becomes in this approximation a system of ordinary differential 
equations. To obtain an accurate description of the problem, the number of classes and thus 
of equations and interaction parameters must be of the order of 10 — 100. Then, the equations 
can be quickly solved with standard numerical software. One is usually interested into the 
large-time behavior of the solutions. It turns out that these are asymptotic equilibrium 
states, which are otherwise virtually impossible to find algebraically, even with powerful 
symbolic software. 

The discretized Boltzmann approach represents a major progress with respect to the 
classical “thermodynamical” method based on a single rational representative agent. We 
could say that in the discretized approach there are many interacting representative agents, 
each one with the average features of a small subset of all individuals. The advantages with 
respect to agent-based simulations are given by a much smaller number of variables and 
by a “portable” analytical formulation independent from the software. The introduction 
of additional heterogeneity, expressed e.g. by further discrete variables beside the income, 
makes the model more realistic, at the price, of course, of an increase in the number of 
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equations. 

If we want to insert a network structure into a kinetic model, we must do so in a proba¬ 
bilistic way, through the concept of link density P(a) (the fraction of nodes having a links) 
and of correlation functions P(j3\a), which give the conditioned probability that an indi¬ 
vidual with a links is connected to one with /3 links. An approach of this kind has been 
proposed by Boguna et al. [151 ]. who have described the diffusion of epidemics through dif¬ 
ferential equations giving the probability for each individual to be infected and in turn to 
infect others. It is well known that in this case the network structure is crucial, and that the 
hubs have a fundamental role in the transmission of the disease. It was found, for instance, 
that on certain networks there exist no minimum infectiousness threshold for the diffusion of 
contagion on the whole network; this is also the case of some viral software in the Interne 


Our kinetic model allows the introduction of a network in a way similar as done in 


m, 


but with some important differences in the structure of the equations, as discussed in Sect. 
El This is in part because in our model for economic exchanges, money is conserved, while 
in a contagion process, unfortunately, the disease can multiply for free. In Sect. IIIII we 
establish some analytical results concerning the existence and uniqueness of solutions of the 
equation system and the conservation of the total wealth. In Sect. [TV]we report the results 
of numerical solutions for the case of scale-free networks and compare them with those of 
previous works on wealth exchange models based on the Bouchaud-Mezard model and on 
agent-based models. Sect. lYl contains our conclusions. 


II. STRUCTURE OF THE MODEL 

The model we discussed in [1] is based on a system of differential equations of the kinetic 
discretized-Boltzmann kind. Society is described as an ensemble of individuals divided into 
a finite number of income classes; the individuals exchange money through binary and 
ternary interactions, leaving the total wealth unchanged. The interactions occur with a 
certain predefined frequency, and several other parameters can also be defined, in order 
to provide a probabilistic representation as realistic as possible. For instance, we can fix 
the probability that in an encounter between two individuals the one who pays is the rich 
or the poor; we can make the exchanged amount depend on the income classes (variable 
saving propensity), etc. After a sufficiently long time the solutions of the equations reach 
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an equilibrium state characterized by an income distribution, which depends on the total 
income and on the interaction parameters, but not on the initial distribution. 

A. Link and income classes, initial conditions, marginal distributions 

As mentioned in the Introduction, we would like to introduce now in our model an ad¬ 
ditional important element of heterogeneity of the population: for this purpose, we suppose 
that different individuals have in general a different number of links of economic nature with 
other individuals. Let this number of links vary between 1 and a maximum value N. The 
population will be divided into N subsets, each one comprising the individuals which have 
the same number of links. This subdivision is regarded as fixed during the evolution, in the 
sense that each individual maintains the same number of links when it becomes richer or 
poorer following the interaction with others. The assumption of constant link distribution 
is clearly only an approximation to a real situation. It can be regarded as adequate for 
situations where the evolution of the income distribution is “soft” in the sense that it does 
not involve massive migrations between classes, which would also likely imply significant 
changes in the connection network. 

When we further subdivide the population into n income classes as usual, the income dis¬ 
tribution of the society at a given instant t is defined by the densities x“(i), with a = 1,..., N 
and i = 1,... , n. The function xf(t) represents the fraction of individuals which belong to 
the i-th income class and have a economic links to other individuals of the population. As 
we will show in Sect. IIII1 the densities xf(t) are normalized at all times: ^2 ai x^(t) = 1. We 
also consider the marginal densities x a [t ) = xf(t) and Xi(t) = ^2 a xf(t). 

The average income of the society (which, due to the normalization of the population, 
also coincides with the total income) is 

N n 

0=1 1=1 

where is the income of the i-th class. As will be shown in Sect. m this is a conserved 
quantity, constant during the time evolution. The average income of the individuals with a 
links is 

n 

n a {t) = ■ 

2—1 


( 2 ) 
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N 

Obviously n = //*(£), but while // is a constant, the a-averages /i Q (f) change in general 

a=l 

in time. Indeed, money circulates among the link classes. It is interesting to compare the 
incomes /z Q (f) at equilibrium, in order to see which link classes a comprise the individuals 
who are, in the average, the richest or the poorest. It is not clear in advance whether a large 
number of links leads, in general, to a larger average income. 

The simplest way to fix the initial conditions is to assume that all individuals are at t — 0 
in the same income class. Previous experience with the model indicates that, for instance, 
if n = 25 and the class incomes grow linearly, then starting off the evolution with all the 
individuals in a class between the 5th and thelOth typically allows a smooth convergence 
towards a unique equilibrium distribution depending only on the total income and exhibiting 
a fat tail. 

Suppose then that xf (0) = 0 for i ^ i 0 , where io is properly chosen as mentioned above. 
According to our previous definitions, we must set x?(0) = P(a), i.e. we distribute the 
individuals at the beginning according to their link classes, all their incomes being initially 
equal. 

More generally, if not all individuals are initially in the same income class, then the 
condition ^Tx“(0) = P{ a ) must be satisfied. 

B. Evolution equations 

In this paper, for simplicity, we shall incorporate a network structure in a particular 
version of our model, namely a version without tax payment and redistribution. The corre¬ 
sponding evolution equations describe only binary exchanges and take the form 

dr(f) n . n 

C hk x h(t)x k (t) ~ C ik x i^) x kif) , (3) 

h,k =1 h,k =1 

where the constant coefficients C l hk , satisfying for any fixed h and k the condition Y^i=i Chk = 
1, express the probability that an individual of the h-th class will belong to the i-th class 
after a direct interaction with an individual of the fc-th class; they define all the features of 

n 

the model, as described in detail in pj. 

We can generalize the evolution equations for the densities xf(t) by introducing the 
network structure as follows: 
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The coefficients C k ° a k ^ G [0,1] in (J4]) express the probability that an individual of the 


h -th class, with a links will belong to the i-th class (maintaining the same number a of 
links) after an interaction with an individual of the fc-th class, with /3 links. They have to 
satisfy 


n 



(5) 


i =1 

for any fixed h, k , a, /3. We take them here as 



( 6 ) 


where the only nonzero elements a)’“ k ^ are those for which h = i, i.e. are 



(7) 


and the only possibly nonzero elements b l h a a k ^ are those of the form 



P((3\a) , 


( 8 ) 


Concerning the terms in ()8j) , 


- we notice first of all that the first line is meaningful only for i < n — 1 and k < n — 1 
and the third line only for i > 2 and k > 2; as for the second line, the first addendum 
on the r.h.s. is effectively present only provided i < n — 1 and k > 2 and the second 


addendum only provided i > 2 and k < n — 1. 

- The coefficients p k ,k £ [0,1] for h,k = 1 express the probability that in an 
encounter between an individual of the /ith-income class and an individual of the fcth- 
one, the one who pays is the /i-individual. As discussed in Q, they carry information, 
in particular, on the level of heterogeneity of individuals belonging to different classes. 
To fix the ideas, we take them here as 


p Kk = mm{r h ,r k }/4r. 
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with the exception of the terms pjj = Vj/2r n for j = 2, ...,n — 1, Ph,i = r\j2r n for 
h = 2 ,..., n, Pn,k = r k/2r n for k = 1 , n - 1 , p hk = 0 for k = 1, n and p hn = 0 for 
h — 1 , n. 


- S denotes the amount of money which is exchanged in an interaction. 


The second-order correlation function P (j3\a) denotes the probability that an indi¬ 
vidual with a links is connected to an individual with p links, ft is normalized as 
YhpP (P\ a ) = 1 and it is a feature of the chosen network structure. It is related to 
the link density P(a) by a “closure condition” 16], which has the form 


aP(P\a)P(a) = pP(a\P)P(P), Vcq P = 1,..., N . (9) 


All this guarantees that the condition (J5]) is verified. 

Note that in the equations used for the epidemic model of 


15] an individual with a 


links is, compared with an individual with one link, a times more likely to contract the 
infection. This is a natural assumption in the case of an epidemic, and highlights the role of 
hubs in the diffusion of the disease. In our model of economic interactions it is possible to 
introduce in (J8]) suitable “amplification” factors, say f a , so that an individual with a links 
has a volume of economic exchanges f a times larger than an individual with one link. This 
could be regarded as a reasonable assumption or not, depending on the example which we 
have in mind and on the context in which we intend to apply the model. If we think of 
the case of the owners of a firm or a store, we conclude that in that case more links mean 
a larger total volume of economic exchanges. In other cases, however, the model without 
amplification is more realistic, implying that an individual with a links has the same total 
economic exchange as an individual with one link, but distributed among the various links. 
This could correspond to the case of a professional who furnishes some service and divides 
his working time among several clients. 


III. SOME ANALYTICAL RESULTS 


We prove here that for the eq.s (j4j) the existence and uniqueness for all t > 0 of the 
solution with prescribed initial conditions hold true. Furthermore, we prove that the total 
wealth remains constant during the evolution. 




Denote 


n N 

£ = {x = {x“}, such that xf > 0 Vi = 1, n, a — 1 ,N and ££*? = i}- (io) 

i =1 o=l 

Theorem 1 Let the elements C(° Q k g (with i, h, k = 1 a,/3 = be as in 

©, ©, ©. For any x(0) G £, the solution x{t) = (x"(f)} of the Cauchy problem for the 
equations ([4j) with initial condition x(0) exists and is unique for all t G [0, +oo). Moreover, 
Vf > 0 one has 

n N 

xf(t) >0 Vf = 1,..., n, Vo = 1,..., N and ££<W = i- (ii) 

i=1 a=l 


Proof: The r.h.s. of the equations © is of class C°° and, in particular, locally Lipschitzian. 
This guarantees the local existence of a unique solution x(t) = {xf(t)} of the Cauchy 
problem. Then, one has xf(t) = xf(0) + f* ds and 

n N n N 

££*?(*) = ££*?(<>) 

i= 1 a=l 2—1 a=1 

r n N n N 1 


+ 




i,h,k=l a,(3=1 


ds 


1 ( 12 ) 


as well, for any t > 0, for which xf(t) exists. This ensures that Yl(=\ Ef=i x ?(t) = 1 as long 
as the solution x(t) exists. 

We denote now 

n N 

G?(x(s)) = J2 C hM x h( s ) x k( s ) ■ 

h,k=l (3=1 

Then, also in view of ©, the equations © can be written as 

= <??(*(*)) - xnt), (13) 

which implies that, if xf(t) is a solution of (fT3lh ■^e t xf(t) = e t Gf(x(t)). Consequently, 

x ?(t) = e _t x“( 0 ) + [ e- t+s G?(x(s)) ds , 

Jo 

and this proves that > 0 for « = 1,..., n and a = 1,..., N, for any t > 0, for which x“(f) 
exists. 
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Since x(t) = {x^(t)} remains within a compact region in R n , its existence for all t > 0 
holds true. | 

We show next that the total wealth expressed by the function 

n N 

W : £ —► R , W{x) = EE r << < 14 ) 

2 = 1 a=l 

is a conserved quantity. Indeed we prove the following result. 

Theorem 2 The value of the total wealth fl 14jf remains constant along any solution of (J3J) . 
Proof: It follows from (J5J), (JSD, dZJ), ((HD and (TTTT) that 


dW 

dt 


n N rlr a 

EE^ 

i= 1 a=l 

n N n N n N n n N 

2=1 a=l k=l (3=1 2=1 a=l h=l k=l 0=1 

n— 1 n— 1 n N N 

h =2 k =2 2=1 a=l (3=1 

n n N N n n N N 


fi X h X k 


n N 

i— 1 a=l 


-YlYlYlY^ r Ka,k,p x i x l + YlYlYlY^ r K^ x n x, k 


a (3 


k =1 2=1 a=l (3=1 
n— 1 n N N 


k= 1 2=1 a=l (3=1 
n—1 n N N 


~x%x? 


+ YlYlY, r K^ X h X 1 + YlYlYjK«,n,fi*h*n 

h=2 i=l a=l/3=1 h =2 i= la=l/3=l 

We will now argue separately on the terms on the r.li.s. in (jT5]h 

• To start with, we notice that, if2</i<n — 1 and 2 < k < n — 1, then 


(15) 


E ii,a ih—l,a , 7/1,a , ih-\-l,a 

r i b h,a,k,p = r h-lb h , a ,k,0 + r hb h ,*,k,p + r ib h+ ia,k,p 

= p k , h S P(a\/3 ) - p h ,k S P{/3\a ) • 


2=1 


Equivalently, the quantity = pk,h S P(a\/3) — ph,k S P(/3\a) is antisymmetric. This 
implies that 

n—1 n—1 N N 

(16) 


n—1 n—1 n N N 

EEEEE’-t 

h =2 fc=2 i=l a=l/3=l 


,*>“ T “ T ^ - 
h,a,k,/3 x hXk ~ 


J2 J2 = 0 > 

/i=2 fc=2 a= 1/3=1 


i.e., the Erst term on the r.li.s. in (fT5j) is zero. 
• Straightforward calculations show that 
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- the second term on the r.h.s. in (fT5j) is equal to 

N N n N N n 

Y^J2( ribl C,k,p x i x l + r ^ a aik ,p x i x l) = 1212 12 sp ( a \P) x i x k ■ ( i7 ) 

a=l/3=lk=l a=l/3=l k= 2 

- the third term on the r.h.s. in (fT5l) is equal to 

N N n N N n~ 1 

Y^J2( rn ~ lb nM X >k + r nK%,k,P< X l) = 1212 12 ~ Pn ’ k S P (P\ a ) X n4 ■ ( 18 ) 

a=l(3=lk=l a=l f3=l k= 1 

- the fourth term on the r.h.s. in (1T5T) is equal to 

n— 1 N N N N n— 1 

EEEK 16 ^^ i+ r hb h c X p x h4) = 121212 -Ph, 1 s P(/3\a)x%x% . (19) 

h= 2 a=l/3=l a=l@=l h= 2 

- the fifth term on the r.h.s. in (JT5j) is equal to 

n— 1 N N N N n— 1 

12 + r h+lbh!a,n,^h4) = 1212 12 p n , h s P(a\/3)x%x p . ( 20 ) 

h=2 £*=1/3=1 0=10=1 h= 2 

Then, adding the four terms in (TT71) — (f2U]l . we get 


a= 1/3=1 


a\p)x“xP — 


1212 PnA S P ( a \P) X l X n = 0 • 

a=l/3=l 


Namely, the sum of the terms from the second one to the fifth one on the r.h.s. in (1151) . is 


zero. 


Summarizing, the r.h.s. in (fl5|) is equal to zero, which proves the claim. 


IV. NUMERICAL RESULTS. DISCUSSION 


We have solved the equations numerically in some simple low-dimensional cases, varying 
the initial conditions and the features of the network. More precisely, we have considered 
scale-free networks with link density P(a) = c/ad, varying the exponent q. We have also 
varied the correlation matrices P({3\a), focusing on the cases of assortative and disassortative 
correlations, but we found that this has little influence on the equilibrium distributions. 

Concerning the dependence on the initial conditions, the results show that the equilibrium 
distributions xf = £“(+ 00 ) all share the following significant properties. 
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1. The marginal income equilibrium distribution x % = ^ a a;“(+ o o) depends only on the 
total income //. 

2. The detailed income equilibrium distribution xf depends on /j and on x a (0), but not 
on the detailed initial condition x"(0). 


Further interesting results concern the difference between the histograms of the quantities 
{xj, i = 1... n}, {xi = 1... n },..., {xj 4 , i — 1... n} which represent the income distribu¬ 
tions of the various “link classes”, i.e. of the classes comprising individuals with 1,2,.. .IV 
links. As mentioned in the Introduction, one could expect that more connected individuals 
generally tend to acquire more wealth. (The strength of this effect is also expected to de¬ 
pend on the possible presence of an “amplified” coupling, compare Sect. ED) It turns out, 
however, that the relative wealth of the link classes is determined by the exponent q of the 
power law of the link density P(ot) = c/a q , the critical value being q — 1. For q < 1 the rich 
individuals tend to concentrate in the classes with many links, for q > 1 in those with few 
links. 

When P(a) — c/a the histograms of the link classes are all similar (Fig. [5]), and more 
exactly one has for any i, j, a, (3 


x. 


X 


x: 


x‘ 


( 21 ) 


or equivalently 


x:i 


X; 


~/3 ~/3 

x i x j 


( 22 ) 


This reminds the energy distribution of a gas containing molecules of different kinds, 
all in equilibrium at the same temperature. In this analogy, individuals with more links 
correspond to molecules with larger cross-section. 

When P(a) = c/a q with q < 1, the histograms of the classes with few links display a 
concentration of individuals in the lowest income classes (Fig. CD). The corresponding Gini 
indices G a (Table DJ are low for small a, because this concentration amounts to smaller 
inequality; in simple terms, in this case the less connected people all tend to be poor. For 
q > 1, on the contrary, it is the most connected people that tend to concentrate in the 
lowest income classes. This appears to mean that when the hubs in the network are very 
few, the mechanism of purely kinetic wealth exchange puts them at a disadvantage. These 
results emerge very clearly from the numerical solutions, although we did not consider a 
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large maximum number N of links yet, since it is not trivial in general to write for any q 
large correlation matrices P(f3\a) which satisfy the network closure condition Q9j) . 

We have mainly considered scale-free networks in our numerical solutions, because there 
is ample empirical evidence that they are dominant in socio-economic interactions. This is 


confirmed also 


interactions 


0 


ly theoretical models in which the network evolves dynamically during the 


. We are not discussing dynamical networks here for two reasons: first, we 
intend to use our model to simulate the effects of fiscal, welfare and macro-economic policies 
on systems which are close to equilibrium, and not in phases of strong growth or economic 
collapse, where the network structure undergoes radical changes; second, we would like to 
compare our results with similar econophysics models, namely the Bouchaud-Mezard model 
and the Yard Sale model, whose network versions also run on a fixed network structure 


13]. In models with dynamical network structure 7-1 111], it is often observed that highly 


connected individuals acquire more wealth, and conversely rich individuals get more and 
more connected. A direct comparison to our case, however, is very difficult, because in 
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11] the rules for the growth of the network are strongly model-dependent, and the rules 


for the economic exchanges contain decisional and game-theoretical elements which can be 
simulated only in an agent-based version. 

Finally note that for fixed link density P(a), the dependence of the behavior described 
above on the assortative or dissortative nature of the matrices P(@_ a) is very weak. 

These results are compatible with those of the cited models [6;, 13j. Bustos-Guajardo and 


Moukarzel |6j] fold (like us) only a weak dependence of the total income distribution on the 
network, but observe a strong dependence on the network of the relaxation time, because 
the network structure has an important effect on the diffusion of wealth in time. We also see 
similar effects in the approach of the numerical solutions to equilibrium (Fig. 2). For q > 1 
the classes with more links tend to lose wealth, starting from equal initial conditions, in favor 
of those with few links. For q < 1 the opposite happens, and for q = 1 there is essentially no 


transfer of wealth between link classes. Loffredo and Garlaschclli m fold that in the case 
of spatially homogenous networks, like those we considered, the total income distribution 
is log-normal; our total distributions are also compatible with a log-normal curve. Loffredo 
and Garlaschclli further observe a transition to a log-normal plus power tail in the presence 
of higher-order correlations, which were not explored in this paper. 
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Link class a 

Gini, q = 1/2 

Gini, q = 1 

Gini, q = 2 

1 

0.24 

0.36 

0.38 

2 

0.33 

0.36 

0.21 

3 

0.38 

0.35 

0.14 

4 

0.40 

0.35 

0.10 

Tot. popul. 

0.37 

0.36 

0.37 


TABLE I: Gini indices of the income histograms of the classes with a links and of the total 
population, in the case of scale-free networks with link density P(at) = c/ofl. 

V. CONCLUSION 


In this work we have set up a consistent mathematical framework, based on kinetic 
differential equations, for the investigation of economic exchanges between individuals linked 
by a fixed network of economic relations. The introduction of a network structure has been 
previously obtained for agent-based models and for the Bouchaud-Mezard stochastic theory; 
it is a novelty in kinetic theory. The use of the network clearly allows a better representation 
of real situations and is an important way of introducing heterogeneity into the model. Our 
results, summarized in Sect. HV]are consistent with some previous results. In addition, our 
formalism allows us to examine the incomes of individuals having a different number of links 
and to test the idea that more connected people tend to become richer. We find that for scale- 
free networks with link density P(ct) = c/a q this is true only if q < 1. A direct comparison 
with models of economic exchanges which include dynamical, disordered or directed networks 
is not possible. Our next objective is not to consider the evolution of the network, but to 
study the effect on a pre-existent network of different economic policies, in line with our 
previous works where the effects of taxation, tax evasion and welfare measures have been 
explored. For that purpose, we will introduce in future work the redistribution terms; we 
decided to omit them here mainly because they would have made the analytical proofs given 
in this paper exceedingly complicated; the main interaction mechanism, however, is already 
contained in the quadratic terms. 
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□ Income, a= 4, q= 1/2 


FIG. 1: Income histograms for link density P(a) = c/a q with q = 1/2. The total population has 
been divided into 10 income classes and 4 link classes. Left: populations of the 10 income classes 
of individuals with 1 link. Right: the same for individuals with 4 links. The fraction of wealthy 
individuals is clearly much larger in the second case. For q = 2, however, the relation is reversed 
(Fig. [3D and for q = 1 the two istograms have exactly the same form (Fig. ED- 





FIG. 2: Time evolution of the total incomes \jP of individuals with a links in the case of link 
density P(a) = c/a q , with q = 1/2, a = 1,2,3,4 (a = 1: upper curve). At the initial time, 
individuals with few links have larger total income because they are more numerous, but then 
they lose wealth in favor of those with more links. Consider, in particular, the individuals with 
1 and 4 links, whose detailed equilibrium income distribution is represented in the histograms of 
Fig. m Since P( 1) = 2P(4), the number of individuals with 1 link is fixed to be twice as large 
as that of individuals with 4 links. At the initial time, their total income is also twice as large, 
because all individuals start in the same income class. At equilibrium, however, the total income 
of the individuals with 1 link has become smaller than the total income of the individuals with 4 
links. This means that the average income of the least connected individuals is now less than half 
the average income of the most connected individuals. More exactly, some of the most connected 
individuals has become very rich (compare Fig. ED, building up a tail in the income distribution, 
while all the individuals with 1 link have got concentrated in the low-income classes. For q = 2, 
exactly the opposite happens (compare Fig. ED- 
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□ Income, a- 4, q=2 


FIG. 3: Income histograms for link density P(a) = c/a q with q = 2. The total population has 
been divided into 10 income classes and 4 link classes. Left: populations of the 10 income classes 
of individuals with 1 link. Right: the same for individuals with 4 links. The fraction of wealthy 
individuals is clearly much larger in the first case. Compare Fig. [1] 





FIG. 4: Time evolution of the total incomes fjf* of individuals with a links in the case of link 
density P(a) = c/a q with q = 2. At the initial time, individuals with few links have larger total 
income because they are more numerous; then those with 1 link gain further wealth from those 
with more links. Compare Fig. [2J 




FIG. 5: Income histograms for link density P(a) = c/a q with q = 1. The total population has 
been divided into 10 income classes and 4 link classes. Left: populations of the 10 income classes 
of individuals with 1 link. Right: the same for individuals with 4 links. The histograms of the 
link classes are all geometrically similar. This reminds the energy distribution of a gas containing 
molecules of different kinds, all in equilibrium at the same temperature. The corresponding Gini 
indices (Table H|) are all equal. 

















































